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Abstract 


We search for three-alpha resonances in by using the complex scaling 
method in a microscopic cluster model. All experimentally known low-lying 
natural-parity states are localized. For the first time we unambiguously show 
that the 0^ state in which plays an important role in stellar nucleosyn¬ 
thesis, is a genuine three-alpha resonance. 
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I. INTRODUCTION 


Carbon, which is the fundamental basis of the chemistry of terrestrial life, is produced in 
red giant stars by burning the helium ash of hydrogen fusion. In order to produce carbon in 
stellar nucleosynthesis, the A = 5 and A = 8 nuclear mass-stability gaps must be bridged. 
Salpeter and Opik pointed out |]l[ that the lifetime of ^Be is long enough, so that the 
a + a ^ ^Be reaction can produce macroscopic amounts of equilibrium ®Be in stars. Then, 
the unstable ®Be can capture an additional a particle to produce stable However, this 
so-called triple-alpha reaction has very low rate because of the low density of ®Be. 

Hoyle argued that in order to explain the measured abundance of carbon in the 
Universe, this reaction must proceed through a hypothetical resonance of thus strongly 
enhancing the cross section. Hoyle suggested that this resonance is a = 0+ state at 

= 0.4 MeV (throughout this paper E,- denotes resonance energy in the center-of-mass 
frame relative to the three-alpha threshold, while T denotes the full width). Subsequent 
experiments indeed found a 0+ resonance in in the predicted energy region. It is the 
second O’*" state (0^) and the second excited state of Its modern parameters Ej. = 0.3796 
MeV and T = 8.5 x 10“® MeV [Q agree well with the old theoretical prediction. 

We mention here that the long lifetime of the ®Be ground state and the existence of 
the 0^ resonance in at the right energy region are only parts of an incredible chain of 
fortunate nuclear coincidences, which makes the abundant existence of carbon and oxygen 
possible. For an interesting account, see Ref. [Q. 

The aim of the present work is to explore the nature of the 0^ state in The facts that 
this state is very close to the three-alpha threshold, and that the alpha particle is strongly 
bound make it probable that the wave function of 0^ has a dominant three-alpha cluster¬ 
ing nature. The low-lying states of including 0^, have been studied in a number of 
macroscopic (with structureless alpha particles) [^,0 and microscopic [§] three-alpha mod¬ 
els. These models reproduce the general features of the low-lying spectrum. However, 
all these models without exception assume three-body bound state- or two-body ®Be -|- a 
scattering state asymptotics for the wave functions. Thus, none of them obeys the physically 
correct three-body boundary condition for states above the three-alpha threshold. Those 
models which use ^Be -f a asymptotics with bound-state-like ®Be [@ are seemingly adequate, 
because the small width of the ®Be ground state makes its wave function very similar to a 
bound state wave function in a large spatial region. However, one must realize that such a 
model predicts the states above the three-alpha threshold, e.g. 0^, to be two-body, ®Be -|- a, 
resonances. This means that currently there is no unambiguous evidence that these res¬ 
onances are intrinsic states of In fact, it was speculated that the 0^ state is not a 
three-body resonance, but an enhancement coming from the > ®Be + a^a + a + a 
sequential decay 0. This idea was, however, criticized by arguing that all the experi¬ 
mental data supported the genuine nature of this state. 

In the present paper we use a method which is able to handle the three-body dynamics 
of the 0^ state correctly. Thus for the hrst time we can unambiguously show whether this 
state is a genuine three-alpha resonance in We also study other low-lying natural-parity 
states of 
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II. MODEL 


Our model is a microscopic three-cluster {a + a + a) resonating group method (ROM) 
approach to the twelve-nucleon system. The trial function of the twelve-body problem has 
the form 


\ 1 > = 


E-4{4”4”4”xS(Pi.P2)}. 

h,h 


( 1 ) 


where A is the intercluster antisymmetrizer, the cluster internal states are translationally 
invariant Os harmonic-oscillator shell-model states with zero total spin, the p vectors are 
the intercluster Jacobi coordinates, li and I 2 are the angular momenta of the two relative 
motions, L is the total orbital angular momentum and [...] denotes angular momentum 
coupling. The total spin and parity of are J = L and vr = (—respectively. 

Putting (|l|) into the twelve-nucleon Schrodinger equation which contains the nucleon- 
nucleon (A^ — N) strong and Coulomb interactions, we get an equation for the intercluster 
relative motion functions y. These functions represent the three-body dynamics of the 
states. In order to determine these functions, we have to use a method which can handle 
three-body resonances. 

There are indirect and direct approaches. The indirect methods, e.g. [^], study the 


three-body problem at real energies, and extract resonance parameters from the three-body 
phase shifts. A similar method has recently been used to study the 0^ state of [0. 
However, the resonant nature of this state and the resonance parameters were extracted 
from a bound state wave function by using the WKB approximation. 

The aim of the direct methods for resonances is to find the complex-energy poles of 
the three-body scattering matrix. For example, in the authors determined the pole 

positions of the three-body R-matrix by analytically continuing the homogeneous Faddeev- 
equation to complex energies. To get a decisive answer to the question of the nature of the 
0^ state, we must use a direct method. 

Our choice is the complex scaling method (CSM) |^. It reduces the problem of asymp¬ 
totically divergent resonant states to that of bound states, and can handle the Coulomb 
interaction without any problem. The main point of the CSM is, that instead of solving the 
original Schrodinger equation for resonances, a new Hamiltonian is defined by 


and the complex equation 


He = U{e)HU-\e), 


He\^e) = 


( 2 ) 


( 3 ) 


is solved. In coordinate space the unbounded similarity transformation U{0) acts on a 
function f{r,f) as 




( 4 ) 


where f describes the angular part of r. For real angles 6, U{0) results in a rotation into 
the complex coordinate plane, whereas for complex 6, U{6) results in a rotation and scaling. 
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Further on, we shall always use real 6 values. In the case of a many-body problem, the 
transformation given by Eq. (|^) has to be performed in each dynamical Jacobi coordinate. 

For a broad class of potentials there is a remarkable connection between the spectra 
of H and Hg [|T^: (i) the bound eigenstates of H are eigenstates of Hg, for any value of 
6 within 0 < 0 < 7 r/ 2 ; (ii) the continuous spectrum of H is rotated by an angle 26] (hi) 
the complex generalized eigenvalues of Hg, e^es = E-^ — iV/2 (with E^, F > 0), belong to 
its proper spectrum, with square-integrable eigenfunctions, provided 29 > |arg£res|- These 
complex eigenvalues coincide with the S'-matrix pole positions. 

In nuclear physics the CSM has been successfully applied to two-body problems, like in a 
RGM description of ®Be [^, in an OCM model of ^°Ne [^] and in the OCM description of 
the resonances of ^°Li |]^. It was also tested for three-body resonances |@, and was used, 
e.g., in a RGM model of ®He, ®Li and ®Be |^, in searching for three-nucleon resonances 




m 


in three-body models of ®He, ^°He and ^^Li and in an RGM model of ®Be and ®B 
Further details and references of the method can be found there. We note, that the 

The latter method was also 


GSM is identical to a contour rotation in momentum space |25 


used, for example, to study three-body resonances in the A = 6 nuclei 

Up to Eq. (^ our treatment of three-body resonances is exact. Since the resonant wave 
functions become square-integrable in the GSM, we can use any bound-state method to 
describe them. We expand the relative motion functions y in Eq. (|^) in terms of products of 
tempered Gaussian functions, exp[-(pi/7j)2]Y)^m^(pi) exp[-(p 2 / 7 i)^]Tz 2 m 2 (P 2 ) (where 

li and I 2 are the angular momenta in the two relative motions, respectively, and the widths 
7 of the Gaussians are the parameters of the expansion), and determine the expansion 
coefficients from the {S'^g\Hg — e|\k 6 ») = 0 projection equation. This way we discretize 
the continuum and select the square-integrable solutions of Eq. (H). We use ten Gaussian 
basis functions in each relative motion. The matrix elements of the complex scaled many- 
body Hamiltonians were calculated in exact analytic forms by using computer algebraic 
techniques. 


III. RESULTS 


In order to avoid any possible model dependence of the conclusions we use three different 
effective N — N interactions. The Minnesota (MN) force was designed to reproduce low- 
energy N — N scattering data [^], while the rather different Volkov 1 (VI) and 2 (V2) 


forces were obtained from htting the bulk properties of s- and p-shell nuclei [^. Each force 
contains an exchange mixture parameter, u and m, respectively. We £x these parameters by 
requiring that the energy of the ®Be ground state be reproduced. The harmonic-oscillator 
size parameters of the alpha particle internal states are chosen to minimize the free-alpha 
energies. Thus, the wave function of the alpha particle is variationally stabilized. The size 
parameters of the alpha particle and the exchange mixture parameters of the N — N forces 
are listed in Table 1 for the three interactions, together with the energies and radii of the 
alpha particle. 

The parameters of the low-lying ®Be resonances, given by the MN, VI and V2 forces, 
are shown in Table 2. For ®Be we use a two-alpha cluster model wave function, similar to 

Eq. 0 
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The relative motion function x is determined by using a two-body scattering approach based 
on the Kohn-Hulthen variational method |^. The a + a scattering phase shifts, coming 
from the MN interaction, are shown in Fig. 1, together with the experimental data. A nice 
agreement is observed. The resulting scattering matrices are continued to complex energies, 
where their poles are localized [p3. The resonance parameters in Table 2 were extracted 


from the complex pole positions. 

This analytic continuation method is in principle equivalent with the CSM, but is nu¬ 
merically far more precise. The CSM has difficulties in localizing very narrow resonances. 
As an illustrative example, we performed CSM calculations for the ^Be resonances and 
found that although the energy of the ground state is well reproduced, the width is strongly 
overestimated. For 6 values which give stable complex energy spectra, with the discretized 
continuum rotated by close to 26, this overestimation is more than two orders of magnitude. 
We could hne-tune the 6 angle and the Gaussian basis parameters to get the experimen¬ 
tal widths, but then the discretized continuum points would be scattered, forming a band 
rather than a line. If we used such a 6 value in the ^^C calculations, the results would be 
disastrous, making the identihcation of the resonances impossible. Therefore, we use such 6 
values in the ^^C calculations which give more or less stable discretized continua. The price 
we have to pay for this is that we cannot resolve very small widths. However, this is not 
a serious problem, because our primary goal is to show the existence of the ^^C states, and 
not to determine the precise resonance parameters. Our model and effective interactions are 
probably inadequate for this latter purpose. 

The N — N interactions, thus set to reproduce the unbound ®Be ground state, are used 
in the ^^C calculations. In this way we ensure that there is no bound two-body subsystem 
in the three-alpha system as shown also by the experimental data. Therefore, our model 
handles the three-body dynamics properly, and the resonances we hnd are genuine three- 
alpha structures. 

We perform calculations for the low-lying natural-parity states of ^^C. In the three-alpha 
wave function ([|) we include /i = 0 for the O’*" and 2+ states, and additionally li = 2 for 
the 1“ and 3“ states. Our test calculations show that the addition of further conhgurations 
hardly influences the results. For instance, if we include the [li,l 2 ]L = [2,2]0 conhguration 
in the 0^ ground state wave function in addition to the [0,0]0 one, we gain less than 1% 
in the three-body binding energy and less than 0.1% in the absolute energy of ^^C. This 
is in sharp contrast to the hndings of macroscopic models which assume structureless 
alpha particles. We believe that the antisymmetrization causes the [0,0]0 state to strongly 
dominate. 

The parameters of the three-body resonances we hnd, using the MN, VI and V2 forces 
are listed in Table 3, together with the experimental values. The ground state of ^^C is 
strongly overbound by MN, roughly reproduced by VI and underbound by V2. Interactions 
VI and V2 give almost the same absolute energy for the ground state, so the difference 
between the energies relative to the three-alpha threshold comes from the different alpha 
particle energies shown in Table 1. However, this is not the case for the differences between 
the MN and VI and V2 ground state energies. The reason of the differences between the MN 
and Volkov results is most probably the different exchange mixture structure of these forces. 
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While the MN force reproduces the deuteron binding energy in an L = 0 space, assumed 
here, the Volkov forces underbind the deuteron. However, the Volkov forces unphysically 
bind the singlet dinucleon states. 

We believe that the ground state should be expected to be overbound in three-alpha 
models. The reason is that the ground state of ®Be is not a perfect two-alpha state. It 
has been shown that the inclusion of ^Li -|- p and ^Be -|- n channels in the ground state 
wave function of ®Be signihcantly increases its “binding energy” |^. Thus, in order to 
reproduce the ®Be ground state energy in a model which contains these conhgurations, the 
N — N interaction (or the a —a interaction in macroscopic models) should be weakened. We 
expect that then the O^" state of would be closer to the experimental position, using the 
MN force, even if the wave function also contains higher-lying rearrangement channels. 
We note, that all macroscopic models underbind the state of p. The only exception 
is Ref. 0 where the “microscopic” a — a potential was constructed in a way which took 
into account some of the effects of the internal structure of the alpha particle. In agreement 
with our MN result. Ref. p found the ground state of to be overbound. 

Table 3 shows that all experimentally known low-lying natural-parity states of are 
reproduced by our model. This means that these states are all genuine three-alpha reso¬ 
nances. In many cases the calculated resonance parameters are far from the experimental 
values. In order to get closer to the experiments, major improvements of the model, e.g., 
the inclusion of rearrangement channels with nonzero spin, would be necessary. 

Our most important result is that we can localize the 0^ state as a three-body resonance. 
In Fig. 2 we show the low-energy part of the O’*" spectrum of the complex scaled Hamiltonian 
with the MN interaction. The two three-body resonances are denoted by the circles. The 
width of the 0^ state is probably overestimated by the GSM (cf. the discussion about the 
GSM and very small widths). The dots in Fig. 2 represent the rotated discretized three- 
body cut. Due to numerical precision problems these points form a band rather than a line. 
We attempted to optimize the Gaussian basis, but could not get a result cleaner than the 
one shown in Fig. 2. The numerical stability of the continuum points is rather sensitive to 
the level of precision maintained during the calculations of the complex many-body matrix 
elements. Although every matrix element was calculated from (very involved) algebraic 
expressions, no special attention was paid to the full optimization of the computational 
algorithm against the loss of numerical precision. Nevertheless, the identification of the 
resonances is unambiguous. 

We found that the parameters of the resonances are approximately independent of the 
rotation angle 6 within a reasonable interval. We performed most of the calculations in 
Table 3 with 6 = 0.2 rad, and checked several of them with 6 = 0.1 rad. Figure 2 shows 
the result of a calculation with 0 = 0.1 rad. We have encountered some problems regarding 
the 3“ states. Using 9 = 0.2 rad, one of these resonances can be lost. That is why we used 
0 = 0.1 rad for this state in Table 3. 

In contrast to the ground state of ^^G, the O 2 state is predicted by all three interactions 
at roughly the same position, and close to the experiment. This is not surprising, because 
the 0^ state, being so close to the three-alpha threshold, is expected to be a more perfect 
three-alpha state than the ground state. 

We would like to note, that although the O 2 state is a three-body resonance, it does 
not decay into an uncorrelated three-alpha final state. Experiments show that the decay 
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proceeds predominantly through the ®Be + a^a + a + a sequential process |^. The 


experimental upper limit for the contribution of the three-alpha decay to the alpha decay 
width of the 0^ state is less then 4%. The dominance of the ®Be -1- a decay over the 3a one 
is the result of the difference between the relative phase spaces 
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IV. CONCLUSION 

In summary, we have studied the resonances of in a microscopic three-alpha model. 
We used the complex scaling method, which allowed us to describe the three-body Coulomb 
dynamics correctly for resonances. We used three different effective nucleon-nucleon interac¬ 
tions and their results are consistent with each other. We have localized all experimentally 
known low-lying natural-parity states in although a better agreement with experiment 
would require major improvements of our model. For the hrst time we were able to unam¬ 
biguously show that the 0^ state of which plays an important role in the astrophysical 
triple-alpha process, is a genuine three-alpha resonance. 
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FIGURES 


FIG. 1. Calculated a + a scattering phase shifts in the center-of-mass frame using the MN 
interaction. Experimental data are taken from Ref. Q. 

FIG. 2. Low-energy eigenvalues of the complex scaled Hamiltonian of the 0+ three-alpha states 
in The dots are the points of the rotated discretized continuum, while the circles are three-alpha 
resonances. The rotation angle is 0.1 rad. 
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TABLES 


TABLE 1. Harmonic-oscillator size parameter (5 of the a internal state, exchange mixture 
parameter u or m of the N — N interaction, and the energy and point nucleon rms radius of 
the free a particle. 



? 

1 

to 

u or m 

Ea (MeV) 

ra (fm) 

MN 

0.6060 

0.93344 

-24.687 

1.36 

VI 

0.5291 

0.57286 

-27.085 

1.46 

V2 

0.5284 

0.60126 

-27.957 

1.46 

Experiment 

— 

— 

-28.269 

1.48 


TABLE II. 
nances in ®Be. 

Energies (relative to the two-alpha threshold) and full widths of low-lying reso- 
All quantities are given in MeV. 


MN 


VI 


V2 

Experiment 


E 

r 

E 

r 

E 

r 

E 

r 

0+ 

0.092 I 

6.15x10“® 

0.092 

2.36x10“® 

0.092 

5.17x10“® 

0.09189 

(6.8±1.7)xl0“® 

2+ 

3.03 

1.39 

2.34 

1.48 

2.26 

1.42 

3.04±0.03 

1.50±0.02 

4+ 

13.10 

4.11 

9.96 

5.89 

9.55 

5.93 

11.4±0.03 

~3.5 


TABLE III. Energies (relative to the three-alpha threshold) and full widths of low-lying natu¬ 
ral-parity three-body resonances in All quantities are given in MeV. 

MN 


VI 


V2 


Experiment 

E 

r 

E 

r 

E 

r 

E 

r 

0+ -10.43’" 


-7.56’" 


-5.27’" 


-7.2746’" 


0.64 

0.014 

0.71 

0.031 

0.83 

0.077 

0.3796±0.0002 

(8.5±1.0)xl0“® 

5.43 

0.92 

4.75 

0.75 

4.68 

0.89 

3.0±0.3 

3.0±0.7 

16.01 

1.74 

15.44 

2.89 

15.91 

3.71 

10.49±0.02 

0.08±0.02 

2+ -7.63’" 


-5.13’" 


-2.47’" 


-2.8357±0.0003’" 


6.39 

1.10 

5.55 

1.11 

5.49 

1.54 

3.89±0.05 

0.43±0.08 

3“ 1.16 

0.025 

1.35 

0.003 

1.85 

0.014 

2.366±0.005 

0.034±0.005 

11.91 

1.69 

12.18 

1.41 

13.80 

2.79 

11.08±0.05 

0.22±0.05 

1“ 3.71 

0.36 

3.72 

0.47 

3.82 

0.72 

3.569±0.016 

0.315±0.025 


^Bound state. 
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